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ABSTRACT 

We explicitly determine the global structure of the SL(2, Z) bundle over the Coulomb 
branch of the moduli space of asymptotically free N = 2 supersymmetric Yang-Mills theories 
with gauge group SU(2) when massless hypermultiplets are present. For each relevant number 
of flavours, we show that there is a curve of marginal stability on the Coulomb branch, diffeo- 
morphic to a circle, across which the BPS spectrum is discontinuous. We determine rigorously 
and completely the BPS spectra inside and outside the curve. In all cases, the spectrum inside 
the curve consists of only those BPS states that are responsible for the singularities of the low 
energy effective action (in addition to the massless abelian gauge multiplet which is always 
present). The predicted decay patterns across the curve of marginal stability are perfectly 
consistent with all quantum numbers carried by the BPS states. As a byproduct, we also show 
that the electric and magnetic quantum numbers of the massless states at the singularities 
proposed by Seiberg and Witten are the only possible ones. 
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1. Introduction and Summary 



Recently, some new methods have been imagined which led to the complete determination 
of the spectrum of BPS states"'" in iV = 2 supersymmetric theories (first in [1] and later in [2]). 
Part of these results had been anticipated at the conjectural level [3,4]. The problem was open 
since the work of Seiberg and Witten [5,6] where it was realized that there may exist some 
regions in the moduli space where the BPS spectrum cannot be continuously related to the 
semi-classical one. This stems from the fact that the duality transformations (monodromies) 
involved in the description of the low energy physics of these theories are not quantum sym- 
metries at all energy scales. In particular, the semi-classical BPS spectrum need not to be 
invariant under the full monodromy group, but only under a subgroup of it, generated by the 
monodromy at infinity M^, which shifts the electric charge of a dyon by an integer multiple 
of its magnetic charge. The discontinuity in the BPS spectrum is possible because there exists 
a curve of marginal stability where usually stable BPS particles become degenerate in mass 
with other BPS states. This phenomena was first highlighted in two dimensional theories in 
[7]. As for the pure Nf = gauge theory, we will show here that in the cases Nf = 1, 2, 3 these 
curves separate what we call a strong-coupling from a weak-coupling region. The latter name 
is used for the region that contains the semi-classical one, where the Higgs vacuum expectation 
value goes to infinity and the theories are weakly coupled. On the other hand, this so-called 
weak-coupling region extends all the way to the curve of marginal stability where the physics 
is actually strongly coupled. The existence of this curve is a genuine non-perturbative effect 
in the theories with zero bare masses. 

In [1] the authors have shown that in the Nf = theory the strong-coupling BPS spectrum 
must be invariant under some duality transformations which do not belong to the monodromy 
group, but which are related to the existence of a global exact quantum symmetry acting on 
the moduli space. There it was a Z2 symmetry coming from the spontaneous breaking of a Zs 
global symmetry down to Z4. Surprisingly, combining this symmetry with the fact that the 
low energy physics cannot be described using a unique set of elementary light fields over the 
whole moduli space (that is, the SL(2, Z) bundle is not trivial), it was shown in [1] that, in the 

\ By BPS spectrum, we always mean the set of BPS quantum states existing in the theory, not the mass 
spectrum. 

§ Note, however, that for large bare masses of the quarks, such curves also exist at weak coupling where 
jumps in the spectrum can be studied using semi-classical methods [8]. 
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strong-coupling region, the BPS states must come in multiplets of the broken symmetry (the "Z2 
pairs"). This was at the basis of the proof that the strong-coupling spectrum consists only of 
the particles responsible for the singularities, namely the magnetic monopole (n e , n m ) = ±(0, 1) 
and the dyon of unit electric charge (in the normalisations of [1,5]), alternatively described 
as ±(±1, 1). Also included in [1] was a new way of determining the weak-coupling spectrum, 
independent of the semi-classical approach. It was shown that no state with a magnetic charge 
\ n m\ > 2 can be present. As noted by Sen [9] this can also be viewed as a consequence of 
the conjectured self-duality of the N = A theory. The same result also follows solely from 
semi-classical reasoning as noted in [10]. The results of [1] have been subsequently confirmed 
in [2] where a completely independent method from string theory was used. 

In the present paper, we will extend and generalize the methods and results of [1] to the 
case where Nf = 1, 2 or 3 hypermultiplets (quarks), without bare masses, in the spin 1/2 
representation of the gauge group SU(2), are included. If not for the bare masses set to zero, 
these are the most general asymptotically free theories of this type. For Nf = 1, the moduli 
space of vacua M.\ is homeomorphic to the compactified complex w-plane, where u = (trc/> 2 ) 
and is the Higgs field in the adjoint representation of SU(2). This is a Coulomb branch on 
which the gauge group is broken down to U(l). For Nf > 2, in addition to a similar Coulomb 
branch, the moduli space has also Higgs branches on which the gauge group is completely 
broken. In this paper we focus entirely on the Coulomb branch where magnetically charged 
particles exist. We will denote this Coulomb branch as M.N f in the following. The expectation 
value a(u) of the scalar field and its dual ajj(u) are given by the two periods of certain 
one-forms on appropriate elliptic curves [6]. As such they satisfy Picard-Fuchs differential 
equations. These Picard-Fuchs equations were derived in [11] where solutions, valid in the 
neighbourhoods of the various singularities, were also given. However, though this was enough 
to determine power series expansions for the prepotential and its dual, no explicit expressions 
valid on the whole Coulomb branch were given. In the following, we will express au and 
a in terms of standard hypergeometric functions and explicitly display the global analytic 
structure of the SX(2,Z) bundle Ejq f over M.N r This will prove useful when discussing the 
strong-coupling spectrum. 

Once we have these solutions, it is indeed straightforward to determine the curves of 
marginal stability where an /a is real. For each Nf = 1, 2, 3, one has a single closed curve 
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diffeomorphic to a circle. Of course, since a BPS state can become massless only at points on 
this curve, as discussed in [1] , the curve goes through all singular points of the Coulomb branch 
(except u = oo). For Nf = 2, the solutions o,d{u), a(u), up to a multiplicative constant, and 
hence the curve are identical to those of Nf = 0. For Nf = 3, the curve is a rescaled and 
shifted version of the Nf = curve. 

To determine the spectra, we will in particular exploit the global Z4_jV/ symmetry acting 
on the Coulomb branch of the moduli space. From this point of view, the richest structure 
is for Nf = 1 which exhibits a Z3 symmetry. For Nf = 3 there is no such symmetry, but 
we will show that there it is not needed for the determination of the spectrum. In this case, 
the strong-coupling spectrum consists of only the dyon* (—1,2) with magnetic charge two 
and the magnetic monopole (0, 1). The weak-coupling spectrum contains the dyons (n, 1) and 
(2n + 1, 2) for all integers n, in addition to the elementary quarks (1, 0) and W bosons (2, 0), 
while we prove that no magnetic charges \n m \ > 3 can exist. For Nf = 2, we have only 
(0, 1) and (1, 1) in the strong-coupling region, the semi-classical spectrum consisting of all the 
dyons (n, 1) where n is any integer, together with the quarks and W bosons. We show that 
no \n m \ > 2 are allowed. For Nf = 1, the strong-coupling spectrum consists again of only 
the three states responsible for the singularities, namely the monopole (0, 1) and the dyons 
(— 1, 1) and (—2, 1). For the weak-coupling spectrum we prove that again no magnetic charges 
\n m \ > 2 can exist, while it contains all dyons (n, 1), n G Z, as well as the elementary quarks 
and W bosons. Finally note that the neutral abelian N = 2 vector multiplet is present in all 
cases over the whole moduli space. 

We begin in Section 2 with a brief overview of the results obtained in [1] for Nf = and 
then discuss some important facts and ideas we will use in this work. There, we also discuss 
in some detail the role of the broken global symmetries and the relevant representations of the 
flavour groups. Then, in Sections 3 to 5, we present our results for respectively Nf = 1,2 and 3. 
These sections are organized as follows: we start by presenting the structure of the singularities 
on the Coulomb branch, and explain why it is unique. The global analytic structure of the 
SX(2,Z) bundle Ej$ f is then determined, and the curve of marginal stability described. 
When it exists, we also study the global quantum symmetry Z4_ jq f acting on M. n } ■ This leads 

* Note that we changed the sign of the electric charge n e with respect to the conventions of Seiberg and 
Witten [5,6]. 
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to the determination of the BPS spectra. We then show that the predicted decay reactions 
across the curve Cjsf f are perfectly consistent with all quantum numbers carried by the BPS 
states, and in particular with the branching rules of the representations of the flavour groups. 



2. Review of Nf = 0, and general considerations 

In this Section, we review or introduce certain results and ideas, which will prove very 

useful in the following. Along the Coulomb branch, the effective action is given in terms of a 

1 d T 

prepotential J- (a), and the knowledge of a(u) and ao(w) = -— — (a(u)) on all of the Coulomb 

— CJUj 

branch A4, which is the (compactified) it-plane, allows one to reconstruct J- {a). Also the mass 
of any BPS state with electric and magnetic charges n e and n m is determined as^ 

m = V2 \an e - aun m \ = V2 \r]((n e , n m ), (ap, a)) \ , (2.1) 

where 77 is the standard symplectic product. In Seiberg and Witten's work [5,6], the phys- 
ical requirement of having certain monodromies around singularities in the moduli space is 
translated by equating dao/du and da/du with period integrals of the only holomorphic one- 
form on an appropriate genus one Riemann surface. The monodromies are determined by the 
asymptotic behaviour of (ajj,a) near the relevant singularities. Near the point at infinity the 
one loop (5 function of the microscopic theory leads to 

a ( N f)( u ) ~ - V2^l 

2 (2.2) 

afff f \u) ~ (4 - Nt) V2^l logw 
An 

and 

Moo = ( Q _ x J • (2.3) 

The singularities at finite u are due to BPS states (n e ,n m ) becoming massless. Using the 
P function of the abelian low energy theory in which (n e ,n m ) is coupled locally, one can 
obtain the corresponding monodromy matrix. When d BPS hypermultiplets (n e ,n m ) become 



f We will equally note the sections of En s by line or column vectors. 
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massless, which occurs when the states carry a d dimensional representation of the flavour 
group, we have 



The easiest way of determining the explicit solution for (op, a) is to use the Picard-Fuchs 
equation satisfied by the period integrals. These have been derived in [11] for the elliptic 
curves given in [6] , and are second order differential equations having two linearly independent 
solutions. The knowledge of the asymptotics fixes the appropriate linear combinations for ajj 
and a. This will be worked out for Nf = 1 and Nf = 3 below. 

We also recall that it follows from (2.1) that for given u the mass of any BPS state is 
proportional to the euclidean length of the vector n e a — n m a]j that lies on the lattice spanned 
by the two complex numbers a(u) and ajjiu). Charge conservation and the triangle inequality 
then imply that a state (n e ,n m ) can only decay if n e and n m are not relatively prime, i.e. 
if (n e ,n m ) = q(n,m) with n,m,q G Z, g ^ ±1. Hence, BPS states with n e ,n m relatively 
prime are stable. This argument fails if ap/a G R, so that the lattice collapses onto a 
single line. On the moduli space M.N f -> the curve of marginal stability C^ f is defined as 



by a path that does not cross the curve C the spectrum of stable BPS states is necessarily the 
same. As already mentioned in the introduction we call the BPS spectra outside and inside 
the curve the weak and strong-coupling spectra, denoted Sw and Ss, and refer to the regions 
as lZ\y an d 71 s- 

2.1. Review of Nf = 

For a pedagogical introduction to this case, see e.g. [12]. We will use the normalization 
convention of [6] for the electric charge, so that n e is always an integer. Indeed, the Dirac 
quantization condition allows states having half the electric charge of the W bosons which now 
have n e = 2. Though this possibility is not realized in the pure gauge theory, for Nf > it 
will be, the elementary quarks (±1, 0) giving an example. Also, a is now defined by (0) = aa 3 . 
Thus with respect to the "old" conventions of [5,1] we multiply n e by 2 and divide a by 2. 




(2.4) 




Hence, at any two points u and u' that can be joined 
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Then (2.1) is still the correct mass formula. The Picard-Fuchs equation is 

^- 1 ^ + i\{ a j (M) = ' (2 5) 

which is equivalent to a hypergeometric differential equation with parameters a — b — — \, 
c = 0. These parameters are not generic but imply that the solutions are of the logarithmic 
type [13] which in turn ensures that one gets SX(2,Z) valued monodromy matrices. The 
requirement that aj)(u) vanishes at u — 1, so that, according to (2.1), the magnetic monopole 
with (n e , n m ) = (0, 1) has vanishing mass at this point, dictates the choice for aj){u). Similarly, 
the asymptotics a{u) ~ y/u/2 as u — > oo determines the choice for a(u). One has 

u — 1/1 1 1 -u\ 

i , . (2-6) 

u + 1\2 „ / 1 1 



2 2 ' «+ 1 

where F is the standard hypergeometric function, see e.g. [13]. Here and in the following, the 
argument of a complex number always runs from — 7r to ir. Thus, ao{u) has a cut on the real 
line from — oo to —1, while a{u) has two cuts, both on the real line, one from — oo to —1 and 
another from —1 to 1, see Fig. 1. 



-1 -1 1 

Fig. 1: The branch cuts of ar>{u) (left) and of a(u) (right). 

The asymptotic behaviour of the functions au(u) and a(u) as u — > oo, 1, — 1 is [1]: 



a D (u) ~ -v / 2w[logM + 31og2-2] 

7T 

a(u) ~ - V2u 



as u — ► oo 



/ \ . u - 1 

a DW 

2 1 u-1, u — 1 1 u — 1 , _ , , 

a(u) ~ log 1 -1 + 4 log 2 

V ; Ti 2vr 2 & 2 2vr 2 V & ; 



as M — ► 1 



G 



aoiu) 



7Y 



a(u) ~ — 



■u + 1, u + 1 u + 1 , tl , 
^-log^- + ^-(l + 41og2) -4 

U + 1, "U + 1 M + l 



log 



+ (— i7r - e(l + 41og2)) +4e 



as m 



where e is the sign of 3m u. The monodromy matrices around the singular points u 
can then be read off from the different asymptotics. One has 



-1 , 

(2.7) 
-1, 1, oo 




Mi 




M_i 




, mU 



3 
-1 



, (2.8) 



where M_i is to be used if the monodromy around u — — 1 is computed with a basepoint in 
the lower half w-plane and M'_ 1 if the basepoint is in the upper half w-plane. This distinction 
comes about since the neighbourhoods of u — — 1 in the upper and lower half planes are 
separated by the branch cuts. This singularity structure corresponds to a massless magnetic 
monopole (0, 1) at u — 1 and to a massless dyon, alternatively described as (2, 1) or (—2, 1), 
at u — —1. 

The curve of marginal stability Cq was studied in [4,14], see Fig. 2. Along this curve, aj^ja 
takes all the values in the interval [—2, 2], with a^/a increasing monotonically from —2 to +2 
as one goes along the curve clockwise from u — —l + ietou = — 1 — ie, with ajjja = at 
u = 0. 



i 


i 

0.86 


< 


> 








-0.86 



Fig. 2: The curve C in the complex w-plane where ^ € R. 
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The mass formula (2.1) implies the following important property: a BPS state with 

— — = r G [—2,2] will become massles at the point u* G Co where — (u*) = r. On the 
n m a 

other hand, the only singularities at finite u on the moduli space are at u = ±1, and hence 
the only BPS states ever getting massless are the magnetic monopole and the dyon. This 
constraint was extensively used in [1] in the derivation of the Nf = spectra. In the Nf > 
case the same constraint will also prove to be very useful. 

Let us finally mention that for Nf = 0, with the present normalisation, the weak-coupling 
(semi-classical) spectrum consists of the dyons (2n, 1), n G Z as well as the W-bosons (2,0). 
The strong-coupling spectrum contains the monopole (0, 1) and the dyon (±2, 1) only [1]. 

2.2. Global symmetries 

Let us recall the field content of the theories we consider. In addition to the N = 2 
vector multiplet containing the gluons, the two gluinos (A, ip) and the Higgs scalar 0, we 
have Nf hypermultiplets (qf , yA<?/?X/) where (qf , qf) are complex scalars and (x , Xf) are 
Weyl spinors. / is the flavour index running from 1 to Nf and (q, x) and (q, x) transform in 
complex conjugate representations of the gauge group. For SU(2) these are actually equivalent 
representations. 

Because N = 2 supersymmetry is not broken in this theory, we have an SU(2)r quantum 
symmetry for which the gluinos and the scalars (q, (p) are doublets. This global exact quantum 
symmetry will play no particular role in the following. At the classical level, since the bare 
masses of the quarks are zero, we also have an U(1)r symmetry under which the gluinos have 
charge 1, the Weyl spinors in the hypermultiplets charge —1, and the Higgs scalar charge 2. 
This symmetry acts on the Higgs scalar and thus potentially on the moduli space. It will be 
of utmost importance in our analysis. Finally, when Nf > 0, which is the case on which we 
focus by now, we have the flavour symmetry which is exceptionally 0{2Nf) here, thanks to the 
pseudo-reality of SU(2) representations. The U(1)r and 0(2Nf) symmetries are anomalous. 
Only a subgroup Z 2 (4_jv / ) of U{1)r survives quantum mechanically. Moreover, the parity 
p in 0(2Nf) is also anomalous, and the quantum flavour symmetry is really SO(2Nf), or 
more precisely Spin(2Nf), as the theory contains spinor multiplets of SO(2Nf) (typically, 
semi-classical reasoning shows that this is the case for the particles with unit magnetic charge 
[6,10,15]). This parity p, which exchanges (q, x) and (q, x) for one flavor of quark and leaves the 
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other fields unchanged, when composed with the U(1)r symmetry, gives rise to an improved 
exact Z>4(4-N f ) quantum symmetry, whose square is nothing but the Z 2 (4_tv/) coming from 
U{1)r alone. 

Z4 symmetry at fixed u 

At a given point u of the Coulomb branch, this 7>4(4-N f ) symmetry is spontaneously broken 
by the vacuum expectation value of the Higgs field down to Z4. This Z4 changes the sign of the 
Higgs field and thus of the electric and magnetic charges. This explicitly shows that when the 
state (n e ,n m ) is present, then the anti-particle (— n e , —n m ) is also present. Thus we will often 
denote these simply by ±(n e ,n m ). However, note that, unlike for Nf = 2, for Nf = 1 and 
Nf = 3 this apparently innocent charge conjugation operation does include a p transformation. 
For Nf = 1, this simply changes the sign of the U(l) flavour charge S. For Nf = 3, p acts 
as the outer automorphism of the corresponding Lie algebra (D3 in the Cartan classification) 
and will change the Spin(6) representation under which the (n e ,n m ) state transforms into its 
complex conjugate. For the spinor representations this amounts to changing the chirality. For 
Nf = 2, particles and anti-particles occur in the same representations of Spin{4). 

Z4_j\r / symmetry on the Coulomb branch 

The spontaneous breakdown of the 7i^_ Nf ^ to the Z4 symmetry is at the origin of the 
Z4_jv / symmetry on the Coulomb branch of the moduli space. This is Z3 for Nf = 1 or Z2 for 
Nf = 2, where the broken generators will relate different points on the Coulomb branch where 
physically equivalent theories must lie. For Nf = 3, no such symmetry exists. The generator of 
this Z3 (for Nf = 1) or Z2 (for Nf = 2) symmetry always contains a p transformation and thus 
in particular will act on the flavour quantum numbers of the states. For Nf = 1 it changes the 
sign of the £7(1) flavour charge. For Nf = 2, the flavor group is Spin(A) = SU(2) x SU(2) and 
p has the effect of interchanging the two SU (2) factor (this amounts to changing the chirality 
for the spinor representations). Note also that in all cases the generator of the symmetry 
acting on the Coulomb branch shifts the 9 angle by 7r and thus should shift the electric charge 
n e of a state (n e ,n m ) by n m . This will be explicitly checked in the following. 

Relations between the quantum numbers 

We list below some constraints on the representation of the flavour group Spin(2Nf) 
a BPS state (n e ,n m ) can carry. These constraints come essentially from the semi-classical 
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approach which allows one to obtain the general form of the wave function for a ijie^ i^m) 
state, valid in the regime where the field theory is weakly coupled. Semi-classical reasoning 
is likely to give the correct answer for the quantum numbers of the states (though not for 
others quantitative features like the physical mass). We will explain in Section 3 how this is 
corrected by non-perturbative effects. For Nf = 1, the flavour group is 5*0(2) = U(l), and 
the abelian charge £ under this U(l) must be such that 2£ = n m — 2n e mod 4. Note that 
this charge £ is not the abelian charge appearing in the central charge Z of the susy algebra 
when the bare masses of the hypermultiplets are non vanishing [6]. For Nf = 2, the flavour 
group is Spin(A) = £77(2) x £77(2) and the irreducible representations can be labeled like 
(2s\ + 1,2^2 + 1) where si and S2 are the spins of the corresponding £77(2) representations. 
Then s\ must be half integer if n e + n m is odd, and must be integer if n e + n m is even. The 
same is true for S2 with n e replacing n e + n m (note that we could interchange the roles of 
si and S2] the convention chosen here and below always correspond to the conventions of 
[6]). For Nf = 3 the flavour group is Spin(6) = SU (4). If 2n e + n m is odd, we have a faithful 
representation of £77(4). The cases 2n e + n m = 1 mod 4 and 2n e + n m = 3 mod 4 correspond to 
complex conjugate representations. We will choose that the monopole (0, 1) is in the defining 
representation 4 of SU (4). If 2n e + n m = 2 mod 4 we have a faithful representation of £0(6), 
and finally if 2n e + multiple of four we have a faithful representation of SU (4) /Z4, 

where here Z4 is the center of £77(4), or the trivial representation. 

We can be more precise for the states having unit magnetic charge. These are always in 
spinorial representations. A state (n e ,+l) has £ = 1/2 or £ = —1/2 (for Nf = 1), it is in 
(2, 1) or in (1, 2) (for Nf = 2), and it is in 4 or in 4 (for Nf = 3), where in each case the first 
entry corresponds to n e even and the second entry to n e odd. Concerning the £ charge of the 
quarks for Nf = 1, note that half of the states (+1, 0) have £ = +1, and the other half have 
£ = -1. 

All these constraints are perfectly consistent with the former analysis of the action of the 
global symmetries. Here we check this only for Nf = 3, and leave the other cases to the reader. 
Note that 2n e + n m = 1 mod 4 is equivalent to — 2n e — n m = 3 mod 4, and thus that charge 
conjugation amounts to complex conjugating the representations in this case, as it should. The 
cases 2n e + n m = mod 4 and 2n e + n m = 2 mod 4 are left invariant by charge conjugation; 
this again is consistent since a faithful representation of £0(6) or £[/(4)/Z4 gives another 
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faithful representation of the corresponding groups by complex conjugation. Let us now check 
the action of the generator of the Z3 symmetry. It indeed acts as charge conjugation as it 
contains a p transformation, since 2(n e + n m ) + n m = —2n e — n m mod 4. 

2.3. Z 4 _jv / SYMMETRY AND BPS STATES 

In this subsection, we will rephrase and sharpen some of the arguments already mentioned 
in Section 4.1 of [1]. The supersymmetry charges, which we will denote by Q a ,l i n chiral 
notation (l</<2foriV = 2 supersymmetry), have charge one under the 7*A{A-N f ) symmetry. 
The relevant piece of the supersymmetry algebra for our purposes is 

{Qa,l,Qp,j} = 2e aP e IJ Z . (2.9) 

This shows that the central charge Z = an e — ar>n m has charge two. Z can be charged under a 
global discrete symmetry group, since such a group is not generated by infinitesimal generators 
belonging to the algebra. The transformation law of Z under the action of the generator of 
the Z^N f symmetry acting on the moduli space is then 

I in 

Z — ► Z' = e 4=7V7 Z _ (2.10) 

Thus we have \Z'\ = \Z\. Moreover, since the Z 4 _jv / symmetry of course maps states of the 
same physical mass into each other {m! — m), we see that a BPS state will be mapped into 
another BPS state: m = y/2\Z\ implies m! = \pl\Z'\. This fact could also be understood in 
terms of conservation of the number of degrees of freedom, as a BPS state lies in a short rep- 
resentation of the supersymmetry algebra unlike the other states with m > \[2\Z\. Moreover, 
if u and u' denote two points related by Z^tv/; and if there exists a matrix G G SX(2,Z) and 
a phase e tuJ determined for instance by (2.10) such that 

Mm^gMvJ^-gM , (2.11) 
\ a J \ a J \n m J \n m J 

the relation \Z'\ = \Z\ will be obviously satisfied. We will show in the following that relation 
(2.11) is indeed realized. 
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To summarize, the existence of a BPS state p = (n e ,n m ) at u implies the existence of 
another BPS state p' = (n' e ,n' m ) at u', u and u' being related by a global symmetry acting 
on the moduli space, with the following properties: they have the same mass, their electric 
and magnetic quantum numbers are related as in (2.11), and the representations of the flavour 
group in which they lie, though they may be different as explained in Section 2.2, have the 
same dimension. 

2.4. A NOTE ON SOME MILD ASSUMPTIONS MADE IN THIS WORK AND SOME OF THEIR 
CONSEQUENCES 

On the Coulomb branch, the coordinate u labels physically inequivalent theories. The 
operators O corresponding to the observables of the theory at u may thus be labeled by u, and 
act in a Hilbert space 7i u . It is natural to suppose that, at least locally, 7i u does not depends 
on u and that the eigenvalues and eigenfunctions of the operators O u vary continuously with 
u. Thus, the mass m(u) of a BPS state, as given in (2.1), must be a continuous function of u as 
one moves in the Coulomb branch. This should also be the case for the electric and magnetic 
charges, but as these are integers, they must be constant. For the same reason, the Witten 
index or its generalizations, used to show that supersymmetry is not broken in the models we 
consider, must be constant. 

This picture cannot be maintained globally. The reason for this is that there exist curves 
of marginal stability, separating strong and weak-coupling regions, allowing the "decays" of 
usually stable BPS states. In [1], it was shown that such decays do happen in the pure 
gauge theory, and one of the most important results of the present paper will be to explicitly 
demonstrate that this is also the case for the 1 < Nf < 3 theories. Thus, at least in the 
BPS sector of the Hilbert space, there must be some kind of discontinuity on this curve. 
The transition between the two sides of the curve may be understood as follows. The theories 
which lie just on the curve should have physically equivalent sectors in their Hilbert space. The 
predicted "decay" patterns on the curve are nothing but the rules for this identification. If this 
is correct, the quantum numbers of the states related by such a decay should be compatible, 
and we will see that this is always the case. Once the identification is made, one can forget 
about the different sectors and only work in a reduced Hilbert space which corresponds to one 
given sector. Now, when we leave the curve, moving into the strong-coupling region, we are 
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left with only the reduced Hilbert space! Typically, the BPS sector of this reduced Hilbert 
space will contain only the states which are responsible for the singularities on the Coulomb 
branch. Note that it is not clear whether the Hilbert space contains other sectors than the 
BPS sector, and if this hypothetic new sectors can undergo discontinuities. 

Another interesting feature of the theories under study is the following. If a spontaneously 
broken global symmetry relates two points u and u' which can be linked by a path which does 
not cross any discontinuity curve, then the broken symmetry predicts the existence of certain 
states in a theory at fixed u. We will show that, in the strong-coupling region, this implies 
that the BPS states come in multiplets of the broken Z^-Nf) symmetry! In the Nf = case 
we had the Z2 pairs of [1], and here we will have Z3 pairs for Nf = 1 and Z2 pairs for Nf = 2. 



3. One flavour of quarks 

3.1. Global analytic structure, and Z3 symmetry 

We begin by discussing the case of a single massless hypermultiplet. We present here all 
the arguments in quite some detail. When dealing with Nf = 2 and Nf = 3 in the following 
sections, we will mainly focus on the new features and try to avoid repetitions. 

As was argued in [6], for Nf = 1 one has three singularities at finite u, related by the 
Z3 symmetry, and due to dyons becoming massless. We will choose the scale and orientation 
on the w-plane so that they are located at the points u\ = e* 71 "/ 3 , 112 = —1 and M3 = e"* 71 "/ 3 . 
Suppose that the singularity at W3 is produced by a state (n e , n m ). The associated monodromy 
matrix is given by (2.4) with d = 1 (here the flavour group is abelian and its irreducible 
representations are all of dimension 1). According to the discussion in Section 2.2, because of 
the Z3 symmetry the states becoming massless at u\ and U2 will be respectively (n e — n rn , n m ) 
and (n e — 2n m ,n m ). A general consistency condition for the monodromy group is actually 



lvl {n e -2n m ,n m ),l M {n e -n m ,n m ),l m {n e ,n m ),l ~ M oo - M(n e -n m ,n m ),l M {n e ,n m ),l M {n e +n m ,n m ),li 

(3.1) 

where Mqq is given by (2.3). This implies n m = ±1. The electric charge is left arbitrary, 
reflecting the possibility of "democracy" transformations, see [1], and we choose n e = 0. This 
shows that the "minimal" choices of [6] are actually the only possible one. 
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Let us now determine the analytic structure of the SX(2,Z) bundle E\. The corresponding 
elliptic curve is [6] 



A 6 

■y = x — ux - 

y 64 



(3.2) 



and da D /du and da^/du are given by its period integrals. Then a D and are given by 

the integrals over appropriate contours 7, of a one-form A, where A is obtained by integrating 
\/2 dx 

with respect to u (up to an exact one-form on the elliptic curve). This yields the 

8n y 

integrals 



2u — 2>x 



^— <£ dx — 

87r I y^ 3 - ux 2 - A§/64 



(3.3) 



for afp and a^. The most convenient choice for the mass scale, leading to singularities on 
\u\ = 1, is such that 



A 



256 
~2T' 



(3.4) 



It is straightforward to show from (3.3), and it was established in [11] that a]-/ and 
satisfy the differential equation* 



(u 3 + 1) 



cr u 
du 2 4 



a (1) 
a D 



(u) = 



(3.5) 



This is not a hypergeometric equation, but it can be transformed into one if we change variables 
[11] to 

v = -u 3 . (3.6) 



Then 



N d 2 / 1 ^ d 1 
v(l - v)—^ + -(1 - v)- 

v J dv 2 3 V J dv 36 



(3.7) 



This is a hypergeometric differential equation with a = b = — i and c 



* It comes from the Picard-Fuchs equation satisfied by the period integrals which are the derivatives of cSp 
and with respect to u. 
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■* 1 

u 2 / 

x 

"3 

3 C 



Fig. 3: Definition of the regions A, B + , B , C and positions of the singularities 
Ul, U2, «3 for Nf = 1 

Note that the mapping (3.6) is 3 to 1 (except at u — 0) and each of the three regions 
A, B (= B + U B~) and C of the -u-plane shown in Fig. 3 are mapped onto the full v -plane. 
The singular points of the differential equation (3.7) are v — 0, 1 and oo, corresponding in the 
■u-plane to 

u Q = , ui = e i7v/3 , u 2 = -1 , u 3 = e~ i7r / 3 , Uoo = oo . (3.8) 

The borders between the regions A, B and C are mapped onto the negative real t>-axis. Any 
solution of eq. (3.7) on the -y-plane will yield a solution of (3.5) in any one of the three regions 
A, B or C. However, to obtain a solution on the whole w-plane we must choose some solution 
in a given region, say in C, and then analytically continue it into the neighbouring regions. 
Note also that although v = is a singular point of (3.7), u = of course is not a singular 
point of (3.5). 

A convenient choice of two linearly independent solutions to the differential equation (3.7) 
are Kummer's solutions U%{v) and Uq(v) [13]: 

Ce(«) = (1 " v)F f|.|,2; 1 - v 
The first solution Uq(v) obviously vanishes as v — > 1, and hence, once appropriately normalized, 
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is a good candidate for a D . The second solution behaves for v — > oo as 

V V6 = (_ u 3)l/6 = ? (3.10) 

where o>(tt) is a region-dependent phase factor: uj{u) = e~ l7r / 2 for -u G 5 + , cj(-u) = e - * 71 "/ 6 for 
tie A, u;(u) = e*^/ 6 for u E C and u;(u) = e 171 "/ 2 for it G Once correctly normalized, C/3 
is thus a good candidate for The asymptotics (3.10) are discontinuous, but as explained 
above, one should take the solution in one region and then analytically continue it to the other 
regions. The correctly normalized function then is 

flWM^^F^-^l;-!) . (3.11) 

This function has one branch cut along the negative real w-axis from the square- root, as 
well as three other cuts extending from the origin of the w-plane to u\, ui and M3 due to the 
hypergeometric function, see Fig. 4. The expression (3.11) for a^\u) obviously has the correct 
asymptotics ~ l^/2u everywhere as u — > 00. 



u 2 



Fig. 4: The cuts of a^(u) 



Next, we determine a D (u). It should be proportional to Uq(v) as given in (3.9) in the 
region A or B or C in which the magnetic monopole becomes massless. Choosing this region is 
a matter of convention, just as for Nf = in [5] it was chosen to have the massless monopole 
at u — 1 (rather than at u — —1). We choose this region to be C with the monopole 
becoming massless at us = e~ l7r / 3 . The normalisation of then is determined from the 
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desired asymptotics (2.2). If we define for all u the function 

/ D ('.) = ^(« 3 +l)f(g|.2;l + 



u 



then for u G C we have 



ag>(«) 



(3.12) 



(3.13) 



Indeed, as u — > oo in region C, the asymptotics of this expression is 



3^ /r— 

— V2m 

47T 



^log(-« 3 ) + ^log2-2 + log3 



3« /r— 

= — V2m 

47T 



4 27T 

log M + -log2-2 + log3 + y 



(3.14) 

in accordance with the asymptotics derived directly from the integral (3.3) [11]. One finds for 
all analytic continuations 



MG B~ : 


a%\u) 


= -/dM + o (1) («) > 


ueC : 


a£\u) 


= e-**f s f D (u) , 


u E A : 


a$(u) 


= e-^ s f D (u)-a^(u 


ueB + : 




= + f D (u)-2a^(u) . 



(3.15) 



We recall that a^^-u) was given in (3.11) for all u. Since the different expressions for a^\u) 
are obtained by analytic continuation through the cuts of F (|, |, 2; 1 + <^w 3 ), which are the 
borders between regions B + and A, A and C, C and B~ , the function a^(u) so defined has 
no cuts at these borders. The only such cut at the border between two regions is between 
B + and B~ . Moreover, in region A, due to the presence of a^\u) in the definition of a^\u), 
there is also a cut from u — to u — u\. The cuts of afp(u) are shown in Fig. 5. 




u 2 



u 3 



Fig. 5: The cuts of a^\u) 
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The solutions (3.15) reflect the Z3 symmetry on the moduli space. One has 

(S ) (e ±2 " /3 «) = (I ) (.) , G W± (3-16) 

provided that u is such that the path t 1— > e ±2z7rf / 3 w, i e [0, 1] does not cross the cut on the 
negative real it- axis, see Fig. 6. 




Fig. 6: The matrices e ±t7T ^Gw± provide the analytic continuations associated 
with the Z3 symmetry. 



We now proceed to give the various asymptotics and monodromy matrices for the singular 
points. First, of course, u = is not a singular point, as we already noticed from the differential 
equation (3.5), although this is not completely manifest on our solution. However, for u — > 0, 
cSp(u) and a^\u) have a power series development of the form cq + c\u + . . ., without any 
logarithms. The constants c$,c\,.. . have different phases on different sides of the cuts going 
through u = 0, as shown in Fig. 4 and 5. Nevertheless, if one takes the asymptotic form 
at some point u close to and analytically continues it along a path surrounding u — one 
obviously gets back the same series, and the monodromy is trivial, as it should. For the other 
asymptotics we find: 



(1)/ x 3i 



4 ZTC 

log M + -log2-2 + log3 + y 



as u 



00 
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a^(u) ~ — — e 

07T 



5^/6 



+ (w-ui) ^5-41og2-21og3- y 
v/2 



> as u — > mi 



07T 



12 + (u — U2) \og{u - U2) 



+ (u- u 2 ) (5 - 4 log 2 - 2 log 3 + me) 



.( 1 )^~-^i±la«(«) + ^c (u-ua) 



> as U — > M2 



(1)( U ) ~ ^ e -W6 

07T 



12e 2i7r/3 + ( M _ U3 ) log ^ _ U3 J 



+ («-« 3 ) (^5-41og2-21og3 + y^ 



V2 2 , 



> as u — > M3 (3.17) 



where e is the sign of ^mu. The monodromy matrices follow as 



-1 3 
-1 



It! 



2 1\ / 1 

-1 •"- = -1 1' 



U2 



1 
-1 2 



1 1V1 U 



'it 2 



3 4 



(3.18) 



where M U2 is to be used for a basepoint in the lower half w-plane (e = —1) and M' U2 for a 
basepoint in the upper half w-plane (e = +1). We have (cf. (3.1)) M Ul M U3 M U2 



Moo = 



M^ 2 M Ul M U:j . As a check one verifies from (3.17) that the particles that become massless are: 
the magnetic monopole (0, 1) at u — us, the dyon (—1, 1) at u — u\ and, at u — 112, the dyon 
described as (—2, 1) if e = +1 or as (1, 1) if e = —1. 
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3.2. The curve of marginal stability 



The curve of marginal stability C\ is defined as the set of all u on the moduli space such 
that w^\u) = a^\u) / a^\u) is real. We already know that this curve has to pass through 
ui,U2 and M3, since massless particles must occur on this curve, and we have w^(u^) = 
0, w^\ui) = —1, etc. As for N* = 0, one can try to determine the curve C\ analytically as 
follows: the Picard-Fuchs differential equation (3.5) implies that satisfies 

{W} = -2-^ (3.19) 

where {/, u} = f'"/f — | (/"//') denotes the Schwarzian derivative of / with respect to u. 
The inverse function u{w^>) then satisfies the differential equation 

u'" 3 (u"\ 2 . ,, 2 u . . 

-v-2 UJ - 2 ^Vtt (3 - 20) 

and it is enough to find the solution for a real argument that satisfies the appropriate 
initial conditions u(0) = u^, u(—l) = u\, etc. 




Fig. 7: The curve of marginal stability for Nf = 1 passes through the three cubic 
roots ui, U2, «3 of —1. It is almost a circle. The numbers —2, —3/2, . . . , 1 indicate 
the values taken by = a D / a^> along this curve. 



But again, the precise determination of the curve C\ is completely irrelevant for our pur- 
pose. Instead, just as we did in [1], we have determined C\ numerically, directly from the 
condition ^mw^iu) = 0. In any case, the Z3 symmetry (3.18) implies that this curve is 
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symmetric under rotations by 4^ around the origin. It turns out to be almost a circle, see 
Fig. 7. 

If we call C~ the portion of C\ between 112 and W3 (see Fig. 8), C° the portion between us 
and u\ and C + the portion between u\ and 112, one can show that 

u G C + G [-2,-1] , u G C° w {1) G [-1,0] , u G C~ <£> w (l) G [0, 1] 

(3.21) 

with tu^^tx) increasing monotonically as u goes around the curve in the clockwise sense starting 
at U2- 

3.3. The weak-coupling spectrum 

Let us first determine the spectrum of BPS states in the weak-coupling region IZw, i-e. 
outside the curve C. 

On the one hand we know that there are the elementary excitations of the perturbative 
spectrum, namely the W bosons (2,0) and the quarks (1,0). Note that it seems that the W 
are degenerate with two quarks and thus may disintegrate as (2, 0) — > 2 x (1, 0). However this 
reaction is impossible here since the W have flavor charge S = while the quarks (1,0) have 
S = 1. Thus the W must be stable in this theory. In the cases Nf = 2 and Nf = 3, such a 
simple argument does not work, and we will then admit the existence of the W in the weak- 
coupling regions, interpreting these states as bound states at threshold, as was done in [6]. 
Clearly, ±(1, 0) and ±(2, 0) are the only stable states of zero magnetic charge. Now, let's turn 
to the magnetically charged states. Since the three states (0, 1), (—1, 1) and (—2, 1) become 
massless at 113,111 and 112, i.e. just on the curve C\, they must exist in both IZw and IZs- The 
monodromy around infinity is a symmetry of the semi-classical weak-coupling spectrum Sw, 
so we conclude that all states ±(n e — 3k, 1), n e — 0, 1, 2, k G Z must be in S\y- But these are 
all the dyons ±(n, 1), n G Z, with unit magnetic charge. It is also easy to see that one cannot 
have dyons with \n m \ > 2 in Sw- If there were such a dyon (n e ,n m ) &Sw then, again, all 
dyons (n e — 3kn m ,n m ) would be in Sw, too. Such a dyon would become massless on C\ and 

lead to a singularity if = 3k G [—2, 1]. There is always such a k G Z. But 

n m n m 

* The distance from the origin is slightly larger at u = ui, ui, 113 and slightly smaller at the borders between 
the regions A, B, C. 
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we know that the only massless dyons have n m = 1. Hence there are no dyons with \n m \ > 2 
in Sw, and 

S w = {±(2, 0) , ±(1, 0) , ±(n, 1) , n G Z} . (3.22) 

Finally, recall from Sect. 2.2 that the states ±(2, 0) have S = 0, half of the quark states ±(1, 0) 
have S — ±1 and the other half have S = =Fl. Concerning the dyons, ±(2n, 1) have S — ±1/2, 
and ±(2n + 1, 1) have S = =f1/2. 




Fig. 8: The definition of the various portions of the curve and of the strong- 
coupling region 



3.4. The strong-coupling spectrum 

Mathematical description of BPS states 

Next, we turn to the strong-coupling spectrum. As for Nf = [1], due to the cuts, see 
Fig. 8, the same BPS state is described by different couples of integers in the different regions 
7Zs+, Us-, T^SO- We call Ss+, &S-, <S,S0 the corresponding spectra of integers (n e ,n m ), 
while Ss means the corresponding spectrum of BPS states (unique throughout IZg). Let's 
first work out these different descriptions. Suppose we have a state described in by 
(n e ,n m ). Now transport it along a path inside the strong-coupling region 1Z$ to a point 
u' ElZs~- Since one does not cross the curve C, the state remains stable and cannot decay. 
The mass m(u) = \f2 \n e a^ l \u) — n m cSp(u) | will vary continuously, even as one crosses the 
cut, since physically nothing happens there. What happens as one crosses the cut is that one is 
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analytically continuing a^\u) and a^\u) to some {u ! ), a^ l \u ! ). The latter are expressed 
in terms of by using the monodromy matrix around W3 as follows: 



ad) 



(i) 



[u 



(3.23) 



Hence 



m(u') = V2 r]({n e ,n m ),{a^\a (1) ){u')) = V2 v{M~ 3 \n e , n m ), {a%\ aV){u')) (3.24) 



and we see that the same BPS state is described in IZs- by 



M, 



u 3 



n r 



n<- 



1T"tn Tie 



in K S - 



(3.25) 



Similarly, in 7Z$+ this state is described by 



71 



n 



M, 



Ul 



n. 



n r 



2n e + n r , 



-n P 



in Ks+ 



(3.26) 



In other words, 



= M^Sso , Sso = M~ 1 S S+ , S s+ = M U1 M U3 S S 



(3.27) 



that is to say 



p = ±(n e ,n m +n e ) £ S S - 



p = ±(2n e +n m , -n e ) G S s + , 

(3.28) 

where p is a unique locally constant section over E\. 

These results deserve some comments. First note that the sign of the couple of integers 
representing p in Ss- and Ss+ is not always fixed* since eq. (3.24) does not fix the sign 
in (3.25) or (3.26). However, in some cases we can determine this sign. For example, the 



* The sign is relevant here since the particle and anti-particle carry opposite flavour charge. 
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monopole state (0, 1) becoming massless at u% can be followed continuously from the weak- 
coupling to the strong-coupling region if one crosses the curve of marginal stability exactly at 
us where (0, 1) is the only charged massless particle and thus is stable. Thus (0, 1) must be 
described by the same couple of integers in Sso and in S$-- The same reasoning applies also to 
the dyon (—1, 1) at u\. Now, the monopole p = (0, 1) need not to be represented again by the 
same electric and magnetic quantum numbers in S$+- Indeed, (3.28) gives p = ±(1, 0) in Ss+- 
This shows that the distinction between electric and magnetic quantum numbers in the strong- 
coupling region is not clear. This is a highly non perturbative phenomenon, and is possible due 
to the non-abelian monodromies. A somewhat similar phenomenon was discussed in [6] when 
considering the deformation of the singularities of the theory with non zero bare mass as the 
latter goes to zero. Nevertheless, the "monopole" (1,0) in Ss+ should not be confused with 
the elementary quark. Actually, (1,0) in S$+ has \S\ = 1/2 while an elementary quark has 
\S\ = 1. The fact that a state (1,0) can have \S\ = 1/2 does not contradict the semi-classical 
constraint 2S = n m — 2n e mod 4, since the strong-coupling region is not continuously related 
to the semi-classical one. The only case where the semi-classical constraints on the quantum 
numbers are also valid in the strong-coupling region is when one can go through the curve C\ 
without any discontinuity as explained above (i.e. in Sso and Ss- for the monopole, but not 
in Ss+)- Actually, when the strong-coupling region needs only be separated in two different 
charts, as is the case for Nf = 2 or 3, this argument implies that the transition functions of 
the SX(2,Z) bundle are compatible with the semi-classical formula for the representations, as 
we will see in Section 4 and 5. But for Nf = 1 we need three charts and at this stage we are 
let with an ambiguity in the signs of (3.28) for general (n e ,n m ). For instance, we do not know 
if (0,-1-1) in Sso an d Ss- is represented by (+1,0) of (—1,0) in S$+, and thus if the flavour 
charge of (+1, 0) in Ss+ is +1/2 of —1/2. We will see in Section 3.5 how to lift this ambiguity. 

Determination of the strong-coupling spectrum 

Now we can determine the strong-coupling spectrum using the Z3 symmetry. We will 
work in Sso an d then obtain Ss± simply from (3.27). So let u ElZso- Recall that by the 
general arguments of section 2.3, there must exist a matrix G\y+ £ SL(2, Z), satisfying (2.11) 
for This matrix was explicitly determined in (3.16). In section 2.3, it was 

shown that the existence of a BPS state (n e ,n m ) at u & 7Zso then implies the existence of a 
state Gw+(n e ,n m ) &Ss+ at u' & lZs+- This same BPS state must also exist at any point 
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in TZso, but there it is described by M~^Gw+{n e , n m ) ESso, according to (3.27). We define 
Gs+ = M~ l 1 Gw+ and conclude that (n e ,n m )ESso implies Gs+(n e , n m )E Sso- In exactly the 
same way one can use the Z3 symmetry with and conclude that there must also 

be a state M U3 Gw-(n e , n m ) ESso- Hence 



Gs+ = M Ui 1 Gw+ = ( „ ) , Gs- = M u . A Gw- = ( _ „ ) = G s ]_ , 



= Gs-^so = Sso > (3.29) 
i.e. iS^o must be invariant under G P S+ for any integer p. Since 

G 3 S± = -1 , (3.30) 

all strong-coupling states come in Z3 triplets: 

G s+ / x G s+ , \ G 5H 

(3.31) 







Remark that these Z3 triplets contain actually six states when one differentiates particles 
and anti-particles. They are multiplet of the broken Z12 symmetry, as announced in Section 
2.4. For example, the Z3 triplet containing the magnetic monopole (0, 1) E Sso also contains 
(1,0) ESso, corresponding to ±(2, — 1) E Ss+ or to ±(1,1) ESs-*, as well as (—1,1) ESso- 
Hence this triplet contains precisely the BPS states that are responsible for the singularities 
and that become massless at 111,112,113 on C\, as one may expect. Now we are going to show 
that this is the only Z3 triplet in Sso- Suppose that (n e , n m )E Sso- If n m = 0, this is either 
(1,0) or (n, 0), n > 2. In the first case it is part of the triplet just discussed. If it were 
(n, 0) then it would be part of a triplet (n, 0), (0, n), (—n,n) in Sso, an d (0, n) would be an 
additional massless state at M3, which we know does not exist. So suppose n m 7^ 0. With 
(n e ,n m ) also (— n m ,n e + n m ) and (— n e — n m ,n e ) must be in S^o- Now one of these three 



★ We explain how to determine the correct signs below. 
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states will become massless on C°, which is the part of the curve C\ in TZso, if one of the three 
expressions 

— = r, -=^ = -_L = ^ l( r), ~ Ue ~ ^ = - f 1 ± -) = Mr) (3.32) 
n m n e + n m 1 + r n e \ r J 

is in the interval [—1,0]. It is easy to see that the functions </?i(r) and (f2(r) are such that 
either r or </?i(r) or (f2(r) is always in [—1,0]. Hence one or the other state of the triplet will 
always become massless somewhere on C\. This contradicts the singularity structure and the 
triplet cannot be in Sso, unless it is the above-mentioned triplet containing (0, 1). Thus 



(3.33) 



Sso 


= {±(0,1) , 


±(- 


-1,1) , 


±(i,o)} 




Ss+ 


= {±(1,0) , 


±(- 


-1,1), 


±(2,-1)} 


= M U1 S S0 


S S - 


= {±(0,1) , 


±(- 


-1,0) , 


±(1,1)} = 





3.5. Disintegrations 

We will now show that all states of Syy can consistently decay into the states of the 
strong-coupling spectrum when crossing the curve of marginal stability. Suppose that the 
curve is crossed somewhere on C°, so that the available states are those of Sso- (The other 
two possibilities can be discussed in exactly the same way and do work out consistently, too). 

First recall the flavour charges S of these states in the strong-coupling region T^o- Since 
the (anti)monopole (0, ±1) in IZso can be obtained by continuous deformation from the semi- 
classical one when crossing the curve C\ at us, it also has S = ±1/2. Similarly, the (anti)dyon 
(=Fl, ±1) can cross C\ at u\ and hence has S = =Fl/2. To determine the charge S of (1, 0) G Sso, 
note that it is described as ±(—2, 1) e Ss+ and as =F(1, 1) G Ss- where it is massless at U2 
and can be connected to the semi classical states. In both cases one obtains S — ±1/2 in Ss+ 
and Ss-- However, as noted above, we do not know whether the state corresponding to (1, 0) 
in Sso is, e -g- i n $S+, (2, —1) with S = — 1/2 or rather (—2, 1) with S = 1/2. Hence we only 
know that (±1,0) has \S\ = 1/2 but we cannot at first sight determine the sign. We will see 
that all decay reactions work out consistently if we assume that (±1, 0) has S = ±1/2, while 
for the opposite choice almost all decays would violate the conservation of the S charge. As we 
know that the decay reactions must take place, as was shown above using arguments involving 
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only the electric and magnetic quantum numbers, we see that we can determine in this way 
the sign of the flavor charge in all cases where an ambiguity remained. Below we list again 
the strong-coupling spectrum, now denoting the states as (n e ,n m )s with the flavour charge 
quoted explicitly. 

S so = {(0,±l) ±1/2 , ( T l,±l) Tl/2 , (±l,0) ±1/2 } 

S s+ = {(=Fl,0)±i /2 , (=f1,±1):fi/2 ,(=F2,±l) ±1/2 } (3.34) 
S s - = {(0,±l) ±1/2 , (±l,0) Tl/2 , ( T l, T l) ±l/2 } . 

In this notation the weak-coupling spectrum contains the states (±2, 0)o, (±1, 0)±i, (±1, 0) T i, 
(±2n, ±l) ±1 / 2 and (±(2n+l), ±l) Tl / 2 , cf. (3.22). All decay reactions across C° must preserve 
the total charges of n e ,n m and S. One indeed has 





, \ /=Fl\ /±1 

< — ► (2n - 1) x + 2n x + An x 

±1/2 V : ' / : 1 2 V : '/ :| 2 V " , 

< — »■ (2n + 2) x [ ° j + (2n + 3) x | 1 

=Fl/2 V :|:1 /=Fl/2 V ±:L /t1/2 



+ (An + 4) x 



and remarkably enough all quantum numbers work out consistently. 



±1 




±1/2 

(3.35) 



It remains to check the mass balance, but this also works out in all cases, just as it did for 
Nf = [1]. For example, for the first disintegration in (3.35) one has on the l.h.s. m = \^2\a^\ 
while the mass on the r.h.s. is (let r denote a^/a^ G [—1,0] since one crosses the curve on 
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C°) 

m = V2 (\a$\ + \a^ ] + a«|) = V2\a^\ (\r\ + |1 + r|) = v^|a^| (-r + 1 + r) = V2\a^\ 

(3.36) 

as it should. 

4. Two flavours of quarks 

4.1. Global analytic structure, Z 2 symmetry and BPS spectrum 

As was argued in [6], we have for Nf = 2 two singularities on the Coulomb branch, 
where BPS particles in a two dimensional representation of Spin(A) = SU (2) x SU (2) become 
massless. Because of the Z2 symmetry, they must be located at two points related by this 
symmetry; we choose the scale of the theory so that the singular points are at u — 1 and 
u = — 1. Suppose that the singularity at u — 1 is produced by a state (n e ,n m ) in the (2, 1) 
of SU{2) x SU(2) (this fixes our convention for the order of the 577(2) factors of 5pm(4)). 
Thus, n e is even and n m is odd, see Section 2.2. The associated monodromy matrix is given 
by (2.4) with d — 2. According to the discussion in Section 2.2, the state becoming massless 
at u — — 1 must then be (n e ±n m , n m ). Thus it lies in the (1, 2) representation, in accordance 
with the fact that Z2 contains the p transformation. The general consistency condition for the 
monodromy group is 

M (n e ,n m ),2 ^(n e +n m ,n m ),2 = M» = M( ne „ nm>nm ) j2 M^^)^, (4.1) 

and this implies n m = ±1. The electric charge is left arbitrary, reflecting the possibility of 
"democracy" transformations, and we choose n e = 0. This shows that, in this case again, the 
"minimal" choice of [6] is actually the only possible one. 

Now, one may note that the monodromies are exactly the same as in the case of the 
Nf = theory in the old conventions of [5] for a and the electric charge. This, together with 
the asymptotics (2.2), completely determines the solution to be 

a^\u)=a(u) 

a D \ u ) = 2 a D(u) , 

where the functions a(u) and an (u) are defined by (2.6). It follows that a^/a^ = ao/2a 
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and that the curve of marginal stability Ci for Nf = 2 is the same as for Nf = (see Fig. 
2). Furthermore, the Z2 symmetry is implemented formally exactly in the same way as in the 
case Nf = discussed in [1]. 

The analysis of the spectrum can thus be done essentially without modification. In the 
weak coupling spectrum we now have the W bosons (2, 0) which are singlets of SU (2) x SU (2) 
and the elementary quarks (1, 0) in the (2, 2) (which is nothing but the defining representation 
of SO (A)), in addition to the dyons (n, 1). The latter are in the (2, 1) or in the (1, 2) according 
to whether n is even or odd. 

In the strong-coupling region, one must introduce two different descriptions of the same 
section p over E2, one for Qmu > (region Ss+) and one for Smii < (region Ss-)- We 
have 

p = (n e , n m ) E S s + P = (n e , n m + 2n e ) G S S - ■ (4.3) 

Note that this transformation law is compatible with the semi-classical formula for the con- 
straints on the representations, since n m + 3n e = n m + n e mod 2. This is related to the fact 
that here the strong-coupling region is separated into only two pieces, as already explained 
in Section 3. One then shows that the states must come in Z2 pairs which are described, for 
example in Ss+, as { ± (n e ,n m ), ±(n e + n m , —2n e — n m )}. Note that the two SU(2) factors 
of the flavour group are interchanged for the two members of a given pair, as it should be. 
Finally, using our by now standard argumentation, one can show that there is only one Z2 
pair in the strong-coupling spectrum containing the monopole ±(0, 1) and the dyon ±(±1, 1). 

5 5+ = {±(l,-l) , ±(0, 1)} = Mi5 5 _ 

(4.4) 

5 5 - = {±(l,l) , ±(0,1)} = M^S S+ . 
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4.2. Disintegrations 



We first list the predicted decay reactions, which here are completely fixed looking only at 
the electric and magnetic quantum numbers. We perform our analysis in Ss+- 




One could now check that these reactions are compatible with mass conservation and spin. 
This is very similar to what was already done in [1]. Here we focus on the flavour quantum 
numbers. Recall that ±(1, —1) is in (1, 2) and that (0, 1) is in (2, 1). The decay of the quark 
is possible since the latter transforms in (2, 2) and we have (1, 2) <g> (2, 1) = (2, 2). The decay 
of the W is also possible since (2, 2) <g> (2, 2) contains the trivial representation. 

What about the dyons? First note that 2® n contains only SU(2) representations of integer 
spin when n is even, and of half-integer spin when n is odd. Thus (2®( n+1 \ 2® n ) will contain 
(2, 1) but not (1, 2) if n is even, and it will contain (1, 2) but not (2, 1) if n is odd. But this 
is exactly what we need for the disintegration of the dyons to be possible, since ±(n, 1) is in 
(2, 1) or in (1, 2) when n is even or odd. 
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5. Three flavours of quarks 



5.1. The structure of the singularities and the global analytic structure 

From the analysis of [6], we know that we have two singularities in the moduli space, one of 
them due to a massless BPS state (n e ,n m ) which is a singlet of the flavour group 577(4) (and 
thus 2n e + n m = mod 4), and the other due to a state (n e , n m ) in the defining representation 
4 of 577(4) (2n e + n m = 1 mod 4). As there is no symmetry acting on the Coulomb branch 
(except CP), the states (n e ,n m ) and (n e ,n m ) are a priori unrelated. Moreover, at the expense 
of shifting u, one can suppose that (n e ,n m ) is massless at u — and (n e ,n m ) is massless at 
u — 1/4 (the scale is then Ag = 64). The consistency condition for the monodromy group is 

M (n e ,n m ),l %i e +n m ,n m ),4 = M » = M (h e ,h m ),4 ^(n e ,n m ),l ■ (5-1) 

This leads to an intricate system of algebraic equations, which have an invariance corresponding 
to the democracy transformations. We did not try to solve this system explicitly. Instead, we 
give an argument which does not require any computations and which fits the line of reasoning 
promoted in this work. 

In the problem at hand, we should again have a curve C3 of marginal stability, diffeomorphic 
to a circle, passing through the singularities (but we do not check the validity of this assertion 
for a general structure of the singularities). The analytic structure will be similar to the other 
cases where only two singularities are present (Nf = 0,2). Because of the monodromy at 
infinity, ajj/a will take all values in the real interval [h e /h m ,h e /h m + 1] = [r, r + 1] along C3. 
Now, we know since the work in [10,15] that all the dyons ±(2n + 1, 2) with magnetic charge 
two do belong to the semi-classical spectrum, and are singlets of £77(4). Since there is always 
an integer n such that (2n + l)/2 e [r, r + 1], we see that this state must become massless 
somewhere on the curve. Of course this state is (n e ,n m ). The choice of an odd integer n e is 
then purely conventional, and we will take (n e , n m ) = ±(—1, 2). Then one may return to (5.1) 
and show that necessarily (n e ,n m ) = ±(—1, 1). This conclude our analysis of the structure of 
the singularities, and shows that the choice proposed in [6] is again the only possible one. 

The elliptic curve is in the case at hand [6] 

y 2 = (x-u)(x 2 -x + u) (5.2) 
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and da^/dw and da( 3 )/dw are given by the corresponding period integrals. With appropriate 
contours 7; on the elliptic curve, cSp and are then given by 



V2 



4 



dx 



2u — x 



8vr J y/ (x — u) (x 2 — x + u) 



(5.3) 



It is straightforward to show [11] that ol^ and given by the integrals (5.3) satisfy the 
differential equation 



z(l 



dz 2 




0, 



(5.4) 



where z can be chosen to be 4u or equally 1 — 4u. This is the hypergeometric differential 
equation with a = b = — ^, and c = 0. We observe that this is exactly the same equation 
as for Nf = 0, except that there z = Hence, we immediately see that two independent 
solutions to (5.4) are given by d]j{8u — 1) and a(8u — 1) with ap and a defined in (2.6). 

It remains to determine the correct normalisations and linear combinations. From the 
required asymptotics (2.2) as u — > 00 together with (2.7) one concludes a^(u) = ^a(8u — 1) 
and 

fK 

(u) = [a D (Su - 1) - 7 a(8« - 1)J . (5.5) 

To determine the constant 7, we simply require to have a massless dyon (—1,2) at u — 1/4, 
that is a(l/4) + 2a D (l/4) = 0. This leads to 7 = 2 and then 



V2 



a$(u) = - — 
16 



cid{su 



V2 
4 



) -2a(8u- 1) 

The positions of the branch cuts are shown in Fig. 9. 



(5.6) 



1/4 

Fig. 9: The branch cuts of a^\u) and a^ 3 \u) 
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The various asymptotics of and at the singularities are 



a$(u) ~ -^-v / 2w[logw + 6 log 2 - 2 + in] 

4"7T 



a (3) (w) ~ \V2u 



as u — > oo 



ag ) ( M )~^[-4 + (4 M -l)log(4 M -l) + (4 M -l)(l-41og2 + m)] { , 

^ f as^-i 

a (3) (w) ~ — [4- (4m- l)log(4u- 1) - (4m- l)(l-41og2)] ' 
Sir 



a^\u) ~ — i(l + e) [-4-4wlog4w + 4w(l + 41og2)] - 

167T 4 

/2 /2 
a (3) (w) ~ — ie[4 + 4wlog4w-4w(l + 41og2)] + — it 

87T 2 



as m — > , (5.7) 



where again e is the sign of 3m u. It is then straightforward to determine the monodromy 
matrices for analytic continuation of (a^\a^) around these three singular points: 

m - ("o -!) • m "«= (J -!) • (-1 !) • M <>= {-I -I) • (58) 

where Mq is to be used if the monodromy around u = is computed with a basepoint in the 
lower half w-plane (e = — 1) and Mq if the basepoint is in the upper half w-plane (e = +1). 

5.2. The curve of marginal stability 

Let's now discuss the curve C of marginal stability. Again this is the set of all u such that 
■u/ 3 ) = Op /a^ 3 ) is real. Since 



,(3) 

^ = % - 4 



a D (8u - 1) _ " 
a(8u - 1) 



= i [«,(8u - 1) - 2] , (5.9) 



this curve is a rescaled and shifted copy of the curve for Nf = 0, which is almost an ellipse, 
see Fig. 10. 
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0.108 








-0.108 





Fig. 10: The curve of marginal stability for Nf = 3 

If we again call C + and C~ the parts of the curve in the upper half and lower half w-plane, 
then 

u> (3) G [-1, -h on C+ , e [-1 0] on C~ (5.10) 

with w^(0 + it) = —1, w^\j) = —\ and w^(0 — it) = 0, as follows immediately from the 
corresponding properties [1] of w we recalled in Sect. 2.1, but also from (5.7). 

5.3. The weak and strong-coupling spectra 

Next, we determine the weak and strong-coupling spectra Sw and Ss, corresponding to 
the regions TZw and 7^5- outside and inside the curve C. In all of 1Z\y, a BPS state can be 
described by a unique couple of integers (n e ,n m ). However, similarly to what happened for 
Nf = 1 or 2, in the strong-coupling region IZs one needs to introduce two descriptions (n e , n m ) 
and (n e ,h m ) for Smti > (in 7Zs+) and Qmu < (in IZs-), for the same BPS state p. By 
exactly the same argument as for the other cases, this time one finds 

p= (n e ,n m ) e S s+ p= (-n e -n m ,4n e + 3n m ) e S S - ■ (5.11) 

Now we are in a position to determine the strong-coupling spectrum. While the semi-classical 
monodromy must be a symmetry of the weak-coupling spectrum, this need not and will 
not be the case for the strong-coupling spectrum. Also, here we do not have any broken 
global quantum symmetry at our disposal. But this is actually not needed. As is by now 
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familiar, a BPS state described by (n e , n m )G Ss+ will become massless somewhere on C + if 
n 1 

— = r G [—1,—-]. This same BPS state is described by (n e ,n m ) &Ss~, and it will be 
n m 2 

massless somewhere on C" if 

ne_ = _ n e + n m = _r±l_ _ 
4n e + 3n m 4r + 3 

is such that (p(r) G [—^,0]. But it is easy to see that the function y?(r) is such that either 
r G (and then tp(r) £ (-5,0)) or (p(r) G [-5,0] (for r ^ (-l,-±)). Hence any BPS 

state in the strong-coupling spectrum becomes massless somewhere on the curve C. But a 
massless BPS state always leads to a singularity. Since there are only two singularities, at u — 
and u — \, the only states in the strong-coupling spectrum are the two states (together with 
their anti-particle) responsible for the singularities, namely the magnetic monopole, described 
as (0, 1) G S S - or (1, -1) G S s+ , and the dyon (-1, 2) G S s ±: 

S S+ = {±(-1,2) , ±(1,-1)} 

(5.13) 

5 S _ = {±(-1,2) , ±(0,1)} . 

It remains to determine the weak-coupling spectrum Sw- As already mentioned, the only 
symmetry here is the monodromy around u — 00, and one must have 

MooSw = S w . (5.14) 

Since the states that become massless at u — and u — \ must be present in the weak and 
strong-coupling spectrum, we know that (—1, 2), (0, 1) G Sw- Applying , k G Z, to these 
states one obtains all dyons (n, 1) and (2n+l, 2), n G Z. The dyons ±(2n+l, 2) obtained in this 
way are all singlets of SU (4) since of course Mqq does not change the flavour representations. 
One might ask whether it is possible to have dyons ±(2n + 1,2)' which are not singlets of 
SU(4) (note that for \n m \ = 1 we are automatically in a spinor representation and thus a similar 
question would be irrelevant). But if this were the case, ±(—1, 2)' would be massless at u — 1/4 
together with ±(—1,2), and thus M\u would be changed, which is excluded. Of course, the 
elementary states, namely the quarks (1,0) and W-bosons (2,0) are also present in Sw- Can 
one have dyons with magnetic charge \n m \ > 3 ? The answer is no. Suppose there were such 
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a dyon (n e , n m )G S\y with n m > 3. Then Sw would contain all dyons (— ) fc M ( ^ ) (n e , n m ) = 

{n e — kn m , n m ). There is always afcGZ such that = k G [—1, 0]. One would 

then conclude that (n e — kn m , n m ) becomes massless somewhere on C. But there are no 
massless states with n m > 3. Hence we conclude that Sw cannot contain states with n m > 3, 
and that 

S w = {±(1, 0), ±(2, 0), ±{n, 1), ±(2n + 1,2), n G Z} . (5.15) 



5.4. Disintegrations 

As in the case Nf = 2, the predicted decay reactions are completely fixed when one takes 
into account the conservation of the electric and magnetic charge. We perform our analysis in 
5s- 




(5.16) 



We focus on the flavour quantum numbers, all other consistency conditions being satisfied, 
similarly to the previous cases. First we present some Clebsch-Gordan series for tensor products 
of the 4 of SU (4) which actually will be enough to obtain the desired results. These are 

4® 2 = 6 © 10 

4® 3 = 20 a ©2 x 20 & ©4 . (5.17) 

4® 4 = 1 © 3 x 15 © 2 x 20 c © 35 © 3 x 45 

We recall that the monopole (0, 1) is in the 4 and its anti-particle (0, —1) is in the complex 
conjugate representation 4, while ±(—1, 2) is a singlet. Now it is clear that the decay reactions 
for the elementary quarks (+1,0) are possible since 4 © 4 contains 6 which is the defining 
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representation of 5*0(6) of which the quarks form a multiplet. Simply taking the complex 
conjugate of the Clebsch-Gordan series show that the corresponding reaction for the anti- 
particles is also possible (this is always the case, and we will thus only study the decays of 
(+2, 0), (n, +1) and (2n+l, +2) in the following). The reaction is possible for the W + because 1 
is contained in 4 <g)4 . Moreover, for n even, 4®( 2n + 1 ) = (4® 4 )®( n / 2 )(g)4 contains 4 as 4® 4 contains 
1, and thus the decay of the dyon of unit magnetic charge and even electric charge is consistent. 
The same is true when the electric charge is odd, because then 4®( 2n +!) = (4®4j®(n-i)/2 ( g | ^3 
and 4 <8>3 contains 4. We end the discussion with the decay of the dyons of magnetic charge 2. 
Here, the desired result follows from the fact that 4® 4 ( n+1 ) contains the trivial representation 
in its Clebsch-Gordan series. 

6. Conclusions 

In this paper, we have shown that the results of [1] for Nf = can be extended to the 
general asymptotically free N = 2 supersymmetric Yang-Mills theories with gauge group SU (2) 
when massless hypermultiplets are present. In all cases, we obtained the global structure of the 
SL(2, Z) bundles, showed the existence and studied the properties of the curves of marginal 
stability, separating the strong from the so-called weak-coupling regions. For Nf = 1 and 
Nf = 2, we used the discrete global Z3 or Z2 symmetry on the Coulomb branch of the moduli 
space and showed that the states in the strong-coupling region must come in multiplets of the 
corresponding spontaneously broken discrete global symmetry Z^_ Nf y We saw for Nf = 3 
that the existence of such a global symmetry was not crucial. Even for Nf = 1 and 2 one 
might not have used it explicitly: it is really the global structure of the moduli space which 
gives all the constraints we need. Of course this global structure is strongly constrained itself 
by the global symmetry when it exists. 

In all cases we determined the strong and weak-coupling spectra of BPS states completely 
and rigorously. The strong-coupling spectra contain precisely only those states that are re- 
sponsible for the singularities at finite points of the Coulomb branch. We showed that the 
discontinuities of the spectra across the curves of marginal stability lead to disintegrations 
that always are perfectly consistent with the conservation of mass, electric, magnetic and 
flavour charges. As a byproduct, we also showed that the electric and magnetic quantum 
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numbers of the massless states at the singularities proposed by Seiberg and Witten are the 
only possible ones. 
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